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(54) Radar level gauge 

* 

(57) The present invention is directed to a method 
and apparatus of operating a FMCW radar tank level 
gauge (14) that measure a level (13) in a tank (10) by 
obtaining a set of data points of mixed transmitted 
waves and received waves. The method comprises 
generating and storing a set of spectral data values 
including & received target marker indicating the level; 
calculating an adaptive set of masking threshold values 
corresponding to at least a portion of the spectral data 
values; comparing the adaptive set of masking thresh- 
old values with the corresponding spectral data values 
to identify at least one spectral data value associated 
with the level; calculating the level of the tank from said 
at least one spectral data value; and reporting the level 
of the tank. 
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Descripti n 

BACKGROUND OF THE INVENTION 

The present invention relates generally to micro- s 
wave radar level gauges. More particularly, the present 
invention relates to radar level gauge circuitry and meth- 
ods of operation used to accurately determine the level 
of a liquid or other material being monitored. 

Radar level gauges are in wide use for making non- 10 
contact measurements of the level of products such as 
process fluids and other materials. These devices utilize 
antennas to transmit electromagnetic waves toward the 
material being monitored and to receive electromag- 
netic waves which are reflected at the surface of the 15 
material being monitored. Typically, a first or reference 
signal having a varying frequency is generated and the 
transmitted electromagnetic waves are produced as a 
function of the frequency of the reference signal. A sec- 
ond signal; is then obtained from the electromagnetic 20 
waves reflected by the surface of the material and 
received by the antenna. The two signals should have 
substantially the same frequency, but different phases. 
A phase shift signal is then generated as a function of 
the phase differences between the reference signal and 25 
the second signal over the range of frequencies. The 
frequency of the phase shift signal is indicative of the 
distance travelled by the electromagnetic waves 
between the antenna and the surface of the material 
being monitored, and thereby of the level of the mate- 30 
rial. Radar gauges of the FMCW (frequency modulated 
continuous wave) type are well known as shown, for 
example in "Mikrowellen messen Fuellstaende M by Prof. 
Dr. Johanngeorg Otto in Design & Elekronik-Sen- 
SOrtechnik 1 3 May 1 997, vol. 1 0/1 997. 35 

Conventionally, radar level gauges utilize complex 
and expensive circuitry to remove, from the phase shift 
signal, frequency components caused by reflect on of 
the electromagnetic waves from surfaces other than the 
material being monitored. For example, the circuitry 40 
attempts to remove interference signals caused by 
reflections from the waveguide, from the sides of the 
tank, and from other surfaces which were not the 
intended target. 

Typically, the interference components cannot be 45 
completely removed from the phase shift signal. Fur- 
ther, when monitoring the level of a liquid experiencing 
turbulent conditions, for example, caused by intentional 
agitation of a liquid, the phase shift signal dill include a 
wide range of frequency components other than the f re- so 
quency components corresponding to the level of the 
liquid which would have been detected in the absence 
of the turbulence. As a result of this and other phenom- 
ena associated with a spectrum analysis of the phase 
shift signal corresponding to a turbulent liquid or to foam 55 
covering the liquid, conventional radar level gauges 
have experienced difficulty in accurately determining 
the level of the liquid in the tank. 



SUMMARY OF THE INVENTION 

It is the object underlying the present invention to 
provide an improved method of operating a FMCW 
radar tank level gauge and a respective gauge that 
measure a level in a tank. This object is achieved by a 
method and apparatus according to the claims. 

BRIEF DES CRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic representation of a tank in 
which a radar level gauge according to the present 
invention is installed; 

FIG. 2 is a block diagram illustrating the radar level 
gauge of the present invention; 
FIG. 3 is a block diagram illustrating in greater detail 
portions of the radar level gauge circuitry of the 
present invention; 

FIG. 4 is a block diagram illustrating various func- 
tions performed and algorithms implemented by the 
radar level gauge circuitry of the present invention; 
FIG. 5 is a graphical representation of a power 
spectral data set. 

FIG. 6 is a state table of operation of a signal proc- 
essor for analysis of power spectral data. 
FIG. 7 is a flow chart of operation of a signal proc- 
essor for analysis of power spectral data corre- 
sponding to the state table. 
FIGS. 8A-8D are plots versus time of the signal-to- 
noise ratio (SNR), the NLF (no line found) value, the 
determined range and the mode of operation, 
respectively, during a series of interfering events. 

DETAILED DESCRIPTION OF THE PREFERRFn 
EMBODIMENTS 

FIG. 1 is a schematic representation of process 
tank 10 in which microwave or radar level gauge 14 of 
the present invention is installed. Process tank 10 is 
filled with a liquid or other material 1 2, the height or level 
of which is to be determined by radar level gauge 14. 
Level gauge 14 is mounted on a tank port 16 at the top 
of the tank, and is sealed relative to the tank. Radar 
level gauge 14 transmits radar energy along a 
waveguide, through port 16 and receives reflected 
energy from the liquid surface 13 to provide an indica- 
tion of the level of the liquid within the tank. Radar level 
gauge 14 is coupled to a remote location 1 1 (for exam- 
ple a control room) via two-wire process control loop 15 
carrying loop current l L . _ 

The current l L in process control loop 15 can be 
representative of the level of the liquid or material in 
tank 10, or of a control signal from remote location 11. 
The current ranges between 4 and 20 mA, according to 
the ISA 4-20 mA standard. Gauge 14 can also be con- 
figured to communicate digitally over the 4-20 mA cur- 
rent, as in the HART® protocol, or can be fully digital 
communications as in Fieldbus. 
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Radar level gauge 14 includes an electronics hous- 
ing 20. a housing to flange adapter 22, a process 
mounting flange 24 and an antenna 26. Radar level 
gauge 14 is mounted to a standoff pipe 28 which is fas- 
tened to the top of tank 1 0, around port 1 6. Tank flange 5 
30 is fastened to standoff pipe 28. Process mounting 
flange 24 is secured with bolts to tank flange 30 and is 
sealed with a suitable gasket. Process mounting flange 
24 supports both adapter 22 and electronics housing 
20. 10 

FIG. 2 is a block diagram illustrating circuitry of the 
radar level gauge in accordance with preferred embodi- 
ments of the present invention. The circuitry located 
within housing 20 of radar level gauge 14 includes com- 
munication circuitry 50, power supply circuitry 55 and rs 
radar level gauge circuitry 100. Communication circuitry 
50 couples to two-wire process control loop 15 for com- 
municating over loop 15. As such, communication cir- 
cuitry 50 transmits information indicative of the level of 
material 12 in tank 10 over loop 15 by controlling the 20 
current on loop 15. In the alternative, as stated above, 
communication circuitry 50 communicates the level 
information on loop 15 digitally over the 4-20 milliamp 
current. Further, communication circuitry 50 can, in pre- 
ferred embodiments, receive control signals from 25 
remote location 11 over loop 15. Communication cir- 
cuitry 50 then provides these control signals to radar 
level gauge circuitry 100 for controlling operation of cir- 
cuitry 100. 

Radar level gauge circuity 100 generates an elec- 30 
trical microwave signal, the reference signal. This signal 
is fed to antenna 26 for transmission into tank 10 by 
waveguides or other microwave transmission lines. 
Antenna 26 receives microwaves reflected off of surface 
13 of material 12 in tank 10. An electrical signal corre- 35 
sponding to the reflected microwaves is provided to 
radar level circuitry 100. Radar level circuitry 100 then 
uses the reference signal and the reflected wave signal 
to determine the level of fluid 12 in tank 10. As will be 
discussed with reference to later figures in greater 40 
detail, radar level gauge circuitry 100 includes numer- 
ous features, some of which are embodied in hardware 
and others of which are embodied in software, which 
allow radar level gauge circuitry 100 to accurately deter- 
mine the level of fluid 12 in tank 10. even under turbulent 45 
conditions. 

Power supply 55 is coupled to each of communica- 
tion circuitry 50 and radar level gauge circuitry 100 and 
provides a source of power, other than power from loop 
1 5, for powering communication circuitry 50 and radar so 
level gauge circuitry 100. However, in some embodi- 
ments, communication circuitry 50 includes its own 
power supply circuitry for generating power from the 
current in loop 15. in these embodiments, power supply 
circuitry 55 can be used to supply power only to radar 55 
level gauge circuitry 100 if desired. 

FIG. 3 illustrates one embodiment of radar level 
gauge circuitry 100 in accordance with preferred 



embodiments of the present invention. As shown in FIG. 
3, circuitry 100 includes signal processor 110. control 
processor 115, quartz oscillator 120, digital-to-analog 
converter (DAC) 125. amplifier 130, low pass filter 135, 
voltage controlled oscillator (VCO) 140. divider 145, fifth 
harmonic generator 150. band pass filter 155, amplifier 
160, amplifier 165, mixer 170, amplifier 175, differentiat- 
ing gain filter 1 80, low pass filter 1 85, amplifier 190, and 
analog-to-digrtal converter (ADC) 195. In some pre- 
ferred embodiments, signal processor 110 is a digital 
signal processor adapted for implementing the various 
algorithms and functions of the present invention. Con- 
trol processor 1 15 is a processor adapted for controlling 
signal processor 110 and for communicating between 
signal processor 110 and communications circuitry 50. 
However, in other embodiments the functions of control 
processor 115 are performed by signal processor 110. 
and therefore control processor 115 can be eliminated. 
In these alternate embodiments, signal processor 110 
communicates directly with communications circuitry 
50, as is shown by the dashed line. 

Quartz oscillator 120 preferably provides a 27 MHz 
reference clock signal for operation of signal processor 
110. Signal processor 110 generates digital signals 
which vary repeatedly over a predetermined range of 
values. DAC 125 converts the digital signals into analog 
signals having corresponding amplitudes. The digital 
signals from system 1 10 are varied to provide a stepped 
voltage signal at the output of DAC 125. The stepped 
analog signals from DAC 125 are amplified by amplifier 
130 and filtered by low pass filter 135 before be pro- 
vided to VCO 140. Filter 135 limits the rate of change of 
the reference signal to limit the bandwidth of the signal 

Sref 

In response to the ramp signal from filter 135, VCO 
140 provides as an output a signal having a frequency 
which varies, preferably as linearly as possible, over a 
desired frequency range. The output of VCO 140 is pro- 
vided both to fifth harmonic generator 150 and divider 
145. Divider 145 converts the VCO frequency into a sig- 
nal in the MHz frequency range. The loop provided by 
DAC 125, amplifier 130, low pass filter 135. VCO 140, 
divider 145 and a counter 143 within signal processor 
110 is used as a feedback loop for linearization of the 
rate of change of the frequency of the VCO output. Fifth 
harmonic generator 150 generates frequency compo- 
nents at multiples of the basic frequency, and is 
designed to be such that the fifth harmonic of the VCO 
signal is particularly salient. Band pass filter 155 selects 
the fifth harmonic, and all other multiples of the VCO 
signal, as well as the fundamental wave, are sup- 
pressed. The level of the signal provided as an output 
by filter 155 is increased by amplifier 160. 

The input to the VCO 140 has a linearly stepped 
voltage as shown at 137. In some preferred embodi- 
ments, the frequency of the reference signal varies 
between a range of about 24GHz to about 26GHz. Ref- 
erence signal S REF is provided, via a waveguide 161, 
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having directional couplers 162 and 163, to antenna 26. 
Antenna 26 transmits the corresponding microwaves 
toward surface 13 of material 12. Antenna 26 also 
receives microwave signals reflected off of surface 13. 
Amplifier 160, also provides reference signal S REF to s 
amplifier 165. Amplifier 165 amplifies the reference sig- 
nal S REF to a magnitude necessary to power mixer 170. 
Amplifier 165 also operates as microwave buffer to pro- 
vide isolation. 

The travel of the electromagnetic waves from 10 
antenna 26 to surface 1 3 and back results in a tank level 
dependent phase difference between reference signal 
Sref and reflected wave signal S RW . Mixer 1 70 receives 
reference signal S REF and reflected wave signal S RW , 
and provides as an output phase shift signal S PS . Phase w 
shift signal S PS is a low frequency signal having a fre- 
quency which is dependent upon the phase difference 
between signals S REF and S RW , and which is thereby 
indicative of the distance travelled by the electromag- 
netic waves and thus of the level of material 1 2. In some so 
preferred embodiments, the frequency of phase shift 
signal S PS ranges from about 330 Hz for a distance trav- 
elled from the antenna to the surface of material 13 of 
about 0.2 meters, to about 50 KHz for a distance trav- 
elled of about 30 meters. 25 

Preferably, mixer 170 is of the type which produces 
a higher frequency output signal as the phase differ- 
ence between input signals S REF and S RW increases. 
Therefore, the further the distance travelled by the elec- 
tromagnetic waves and thus the lower the level of mate- 30 
rial 13, the higher the frequency of phase shift signal 
S PS . Since the amplitude of the electromagnetic waves 
decreases as the distance travelled increases, the 
amplitudes of signals S RW and S PS also tend to 
decrease as the level of material 13 decreases. Con- 35 
ventionally, this often required complex AGC (automatic 
gain control) circuity to compensate for the decreased 
amplitude of Higher frequency phase shift signals. Fur- 
ther, when monitoring the level of turbulent liquids, the 
AGC circuit complexity necessarily increases. In addi- 40 
tion to increasing the cost of the level gauge, the com- 
plex AGC circuitry decreases the reliability of the 
electronics in conventional level gauges. 

The amplifier 175 is coupled to the output of mixer 
1 70. Amplifier 1 75 receives phase shift signal S PS and <s 
amplifies it to match the level for optimum use of the 
analog-to-digital converter 195. 

The input of differentiating gain filter 1 80 is coupled 
to the output of amplifier 175 such that differentiating 
gain filter 180 receives phase shift signal S PS from- the so 
amplifier 1 75. Differentiating gain filter 1 80 applies a fre- 
quency dependent gain to phase shift signal S PS and 
provides the amplified phase shift signal at the output of 
filter 180. In preferred embodiments, as the frequency 
of phase shift signal S PS increases, the gain applied to ss 
the phase shift signal increases as well. Thus, higher 
frequency phase shift signals, which have lower ampli- 
tudes as a result of the amplitude loss of the electro- 



magnetic waves as they travel further to and from 
material 13, are amplified more than are lower fre- 
quency signals. Also, preferably the gain of differentiat- 
ing gain filter 180 varies linearly as a function of the 
frequency of the phase shift signal over a desired fre- 
quency range, for example over the range from about 
100 Hz to about 50 KHz. 

Inclusion of differentiating gain filter 180 in accord- 
ance with preferred embodiments of the present inven- 
tion eliminates the need for costly AGC circuit 
components. By applying a frequency dependent gain 
to the phase shift signal, the need for other AGC is elim- 
inated or minimized. Further, differentiating gain filter 
180 provides improved responsiveness as compared to 
many AGC circuit designs. This is very advantageous 
when monitoring th£ level of turbulent liquids. Also, 
unlike AGC circuitry which may need to be calibrated for 
the particular material being monitored, differentiating 
gain filter 180 can be used with virtually any material 
without recalibration. Differentiating gain filter 180 is 
known in the electronics art, for example from Horowitz 
et al., The Art of Electronics, Second Edition, Section 
4.20 (1989). 

Phase shift signal S PS is next provided at the input 
to low pass filter 185. In preferred embodiments, low 
pass filter 185 has a cutoff frequency above the mean- 
ingful frequency range of phase shift signal S PS , for 
example approximately 50 KHz, in order to fflter out 
noise. The output of low pass filter 185 is coupled to 
amplifier 190 for amplification of signal S PS . Finally, 
ADC 195 digitizes phase shift signal S PS and provides 
the digital signal (set of data points of mixed transmitted 
and received waves at each transmitted frequency) to 
signal processor 1 10 for storage. Signal processing 110 
is preferably a microprocessor based circuit adapted to 
receive signal S^ and provide as an output to commu- 
nications circuitry 50 a signal or information indicative of 
the level of material 13. The functions and algorithms 
implemented by signal processor 110, some of which 
can be embodied in hardware and some of which can 
be embodied in software, are discussed below in 
greater detail. 

FIG. 4 is a block diagram illustrating functional 
steps performed by signal processor 110 in preferred 
embodiments of the present invention. As illustrated in 
the functional blocks or steps of FIG. 4, signal processor 
1 10 receives the digitized phase shift signal from ADC 
195 and stores the digitized phase shift signals in ran- 
dom access memory (RAM) as indicated at step 210. 
Next, at step 240, the stored phase shift signal data is 
analyzed. The various steps and features of the analy- 
sis of the signal data illustrated in step 240 are 
described below in greater detail with reference to 
FIGS. 5-8D. After analysis of the signal data at step 240, 
signal processor 110 provides to communication cir- 
cuitry 50 information indicative of the level of the mate- 
rial in tank 10. As described previously, communication 
circuitry 50 reports this level indicative information over 
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two-wire loop 15, or to a local display. 

FIG. 5 illustrates pictorially one set of power spec- 
tral data 300 obtained from FFT (fast fourier transform) 
in step 240 in FIG. 4. Preferably, before FFT transforma- 
tion, the time signal data stored at step 21 0 is Hamming 
weighted. Generally, the power spectral data 300 com- 
prises an integer number of frequency components 302, 
wherein the magnitude of each frequency component 
302 is obtained and stored in suitable memory "bins". In 
a preferred embodiment, the number of bins equals 
51 2, which corresponds to frequency components in the 
range from 0 Hz to approximately 50 KHz. Generally, 
the signal processor 110 processes the power spectral 
data in order to locate a target set of frequency compo- 
nents 304 that is indicative of the level of liquid surface 
13 in the process tank 10. In view that the radar level 
gauge 14 is commonly mounted at the top of the proc- 
ess tank 10, the level of liquid surface 13 in the process 
tank 10 can be expressed as a function of the range of 
the liquid surface 13 from the antenna 26 of the radar 
level gauge 14 since the location of end of the antenna 
26 from the bottom of the process tank 1 0 will be known. 

Briefly, the power spectral data 300 is processed in 
order to locate a frequency component 305 with the 
largest magnitude (target marker) in a greater extent of 
he frequency spectrum of the power spectral data indi- 
cated by double arrow 308. Frequency components 302 
adjacent the largest frequency component 305 define a 
lesser extent of the power spectral data or a "window" 
312, which is used for further, processing. The limits of 
the window 312 are indicated at w. and w + in FIG. 5. 
Using the window 312 allows for efficient signal 
processing of subsequent sets of the power spectral 
data wherein the lesser extent of the frequency spec- 
trum for each successive frequency spectrum is proc- 
essed in order to calculate the level of liquid surface 13 
in the process tank 10. 

To locate the frequency component 305 having the 
largest magnitude in the power spectral data 300, an 
"interference buffer 320 is calculated. Preferably, the 
interference buffer 320 is calculated continuously for 
each new set of power spectral data 300 obtained, if a 
maximum frequency component is found. The interfer- 
ence buffer 320 comprises a value for each frequency 
component 302 or bin. As explained below, the value 
calculated for each frequency component 302 can be a 
function of any one or, preferably, a combination of the 
following: the current magnitude of the frequency com- 
ponent 302, the maximum magnitude recorded for the 
frequency component 302, (if desired, within a selected 
time period), other previous values of the frequency 
component 302 within a selected time period, and a 
function of values of adjacent frequency components 
302. A set of masking threshold values 324 is a type of 
interference buffer where the values have been 
increased by a fixed constant, for example 6dB. 

In a preferred embodiment, values for the interfer- 
ence buffer 320 outside of the search window 312 can 
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be adaptively calculated as follows: 

1) by adapting individual interference buffer values 
as a function of the associated spectral data value 
and at least one adjacent spectral data value. This 
can be expressed as: 
ave[n]=(1/(2*avejen+1 - 

n_exdude))E i . nMWB _ to „ ,H « we - ,en maonIH 
for i g {tpos-twin...tpos+twin} 

and 

m[n]=MAX(ave[n],magn[n]) 
where 



n 

avejen 

n_exclude 

magn 
tpos 

twin 



bin number in the interference buffer 
is one-half of the width of the averag- 
ing window (bins) 

= #({n-avejen.. .n+avejen}n {tpos- 
twin...tpos+twin)) 
values of power spectral data 
bin number having the largest magni- 
tude in the range r. to r + 
is target window, width = 2*twin+1 
(bins) 



and/or, 2) by adapting individual interference buffer 
values as a function of the associated spectral data 
value and past interference values. This can be 
expressed as: 

ib[n,t]=m[n,t](1/M)+ib[n f t-1](M-1)/M 

ib interference buffer 

t time 

M time constant 

Using the aforementioned criteria, the interference 
buffer 320 is a set of values which when compared with 
the current power spectral data provides less drastic 
changes in values than the current set of power spectral 
data. This is illustrated by the dashed line 320 in FIG. 5. 
It should be noted that values of the interference buffer 
320 within the window 312 are not based on the current 
values of the frequency component 305 and adjacent 
frequency components 302 within the window 312, but 
rather, can be a simple linear interpolation between the 
values of the interference buffer 320 at the limits of the 
window 312 (w. and w + ). 

Once the values of the interference buffer 320 have 
been determined, a set of masking threshold values 
indicated at 324 is calculated. In the embodiment illus- 
trated, the set of masking threshold values is a simple 
scaler quantity, for example 6dB above each value in the 
interference buffer 320. The set of masking threshold 
values 324 provides a filtering function so that slight 
changes in values of successive sets of power spectral 
data 300 do not get viewed as a possible target. In con- 
trast, valid target sets, for example, target set 304 in the 
window 312 has values which exceed corresponding 
values of the set of masking threshold values so that 



5 



BNSOOCID: <EP 08876S8A1. 1 > 



9 



EP 0 887 658 A1 



10 



these components of the power spectral data 300 can 
be properly processed in order to determine the level of 
the liquid surface 13 in the process tank 10. 

In one preferred embodiment, analysis of the power 
spectral data 300 is performed using two general 
modes of operation. Referring to FIGS. 6 and 7, a 
"search mode" is indicated generally at 340. The search 
mode 340 includes process steps where the signal 
processor 110 analyzes the greater extent of the fre- 
quency spectrum 308 looking for a target set 304 above 
the sea of masking threshold values 324. Upon location 
of a confirmed target set, the signal processor 110 
enters a "track mode" generally indicated at 342. During 
the tracking mode, the signal processor 1 10 locks onto 
the target set and creates the window 312 for purposes 
of analysis in successive sets of the power spectral data 
300. The signal processor 110 will remain in the track 
mode 342 until target sets are determined to be lost 
within the window 312. If no target sets are found in the 
window 312, the signal processor 110 will return to the 
search mode 340. Thus, if line_track_mode in Fig. 6 
equals TRUE, this implies "track mode" whereas, if 
line_track_mode equals FALSE, this implies "search 
mode". 

Initially, the signal processor 110 enters the search 
mode 340 with an initial set of masking threshold values 
324 as indicated at block 350. Preferably the initial set 
of masking threshold values includes values, which are 
associated with static reflections from the antenna 26 or 
other fixtures associated with the top of the process 
tank 10. It is assumed that the level of liquid surface 13 
in the process tank 10 will not exceed level r_. In a pre- 
ferred embodiment, r. is adjustable by the user. 

It is assumed that the level of the liquid surface 1 3 
of the process tank 10 will not be measured within the 
range to the right of r + . The value r + is also preferably 
adjustable by the user and is used to fitter out fixtures at 
the bottom of fhe process tank 10 such as agitators or 
the like, or when the bottom of the process tank 10 is 
curved. All values of the initial set of masking threshold 
values within the range r. to r + are set at a suitable value 
so that frequency components 302 within this range are 
indicative of possible targets. 

After creation of the set of initial masking threshold 
values, program flow continues to block 352 where a set 
of power spectral data 300 or spectrum similar to FIG. 5 
is obtained. As discussed above, the signal processor 
110 then analyzes the power spectral data 300 in the 
range of r. to r + to locate the frequency component 305 
with the largest magnitude that exceeds the corre- 
sponding value of the initial set of masking threshold 
values. If such a frequency component is not found as 
indicated at block 354, program flow returns to block 
352 where another set of power spectral data 300 is 
obtained. If such a frequency component is not found 
within a predetermined number of successive iterations, 
the signal processor 1 1 0 will enter an alarm state where 
a suitable alarm signal is provided through the commu- 



nication circuitry 50 to the remote location 1 1 (FIG. 1). 

If a useable frequency component 305 is found as 
indicated at block 354, program flow proceeds to block 
360 where the interference buffer 320 is updated. As 

s stated above, the interference buffer 320 is used to cre- 
ate a new set of masking threshold values 324. A coun- 
ter "LF" at block 362 is used to count the number of, 
preferably, successive occurrences of a suitable fre- 
quency component in successive sets of the power 

to spectral data. When the counter LF has reached a pre- 
determined value at block 364, for example eight (which 
can be adjustable in the field), program flow enters the 
track mode 342. 

In a preferred embodiment, when a target set is 

75 found at block 354 a window 312 is subsequently cre- 
ated at block 360, the window 312 is then used succes- 
sively until the counter LF exceeds the predetermined 
value. If a target set is not found in the window 312 
before the counter LF has exceeded the predetermined 

20 value, the counter LF is reset and again the next search 
will be over the range r. to r + . 

Upon entry into track mode 342, the signal proces- 
sor 110 uses Discrete Fourier Transformation (DFT) 
within the window 312 to obtain a range value indicative 

25 of the level of liquid surface 13 in the process tank 10. 
The calculated level of the liquid surface 13 in the proc- 
ess tank 10 is provided to the communication circuitry 
50 in FIG. 3 to be transmitted to the remote location 1 1 
(FIG.1). 

30 Preferably, program flow will remain in the track 
mode 342 indefinitely as long as a set of targets of suf- 
ficient magnitude is located in each successive set of 
power spectral data 300 or spectrum. After calculation 
of the output value using DFT at block 370, program 

35 flow continues to block 372 where another set of power 
spectral data 300 is obtained. Analysis is performed at 
block 374 on the power spectral data 300 within the win- 
dow 31 2 (w. to w + ). If a target set is found within the win- 
dow 312, program flow continues to block 376 where an 

40 interference buffer is created and a new set of masking 
threshold values are calculated similar to block 360 
described above. Program flow then continues back to 
block 370 where DFT analysis is used to determine a 
new output value indicative of the level of the liquid sur- 

45 face 13 in the process tank 10. In a preferred embodi- 
ment, computation of the output value indicative of the 
level 13 includes a coarse computation followed by a 
fine computation to provide a more accurate value. The 
fine computation can be performed using DFT or Chirp 

so z-Transformation. 

If a target is not found in the window 31 2 for a set of 
power spectral data at block 374, program fbw contin- 
ues to block 382. A second counter "NLF" is used to 
count the number of occurrences which a target set was 

55 not found in successive sets of power spectral data at 
block 382. Block 384 determines if the counter NLF has 
exceeded a predetermined value, for example 20 (which 
can be adjustable in the field). Assuming at this point 
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that the counter NLF has not exceeded this predeter- 
mined value, program flow preferably continues to block 
382 where the counter NLF is incremented by one. 
Another set oi power spectral data is then obtained at 
block 372. In a preferred embodiment as illustrated, 
before a new set of power spectral data is obtained at 
block 372, the window size increases in order to search 
a slightly larger extent of the frequency spectrum. In a 
preferred embodiment, the rate in which the window 
increases in size, is given by the following: 
The size of the search window used depends on the 
value of nlf counter: 

swin = max(lf_window_minsize, roc_bins * max(0,nlf 

counter-delay)) 

where 

lf_window_minsize 

the minimum size of the window, if target found 
roc_bins=roc/delta_r(delta_r 

is the rangejgate in meter (distance between two 

bins)). 

roc 

rate of change 
delay 

delay before the window is enlarged 

Thus the window increases over time with the rate 
given by the rate-of -change roc_bins when no target is 
found. If a target set is located in the window 312, the 
counter NLF can be reset. f 

If the counter NLF exceeds the predetermined 
value, program flow continues to block 390 where anal- 
ysis is performed to determine whether the process 
tank 10 is empty or there has been a loss of signal. 

FIGS. 8A-8D illustrate an example of the manner in 
which signal processor 110 of radar level gauge 14 
functions in its two modes of operation, the search 
mode and the tracking mode, as described above with 
reference to FIGS. 6 and 7. The example illustrated in 
FIGS. 8A-8D corresponds to operation of radar level 
gauge 14 measuring level in a tank partially filled with 
peanut oil, and under several adverse conditions. To 
best illustrate the manner in which signal processor 110 
functions in these two modes of operation, there are 
provided in FIGS. 8A-8D four plots which illustrate, 
respectively, signal-to-noise ratio (SNR) FIG. 8A, the 
value of the no line found (NLF) variable FIG. 8B, the 
actual range or distance from the top of the tank to sur- 
face 13 of material 12 as measured with an independent 
sensor FIG. 8C, and the mode (search or track) of oper- 
ation, each plotted against the same time scale FIG. 8D. 
The time scale of each of FIGS. 8A-D are represented 
by numbers of data points, with each data point repre- 
senting a time interval between six and seven seconds. 

Shortly after the initial time period, the signal proc- 
essor enters the tracking mode as illustrated in FIG. 8D. 
As illustrated in FIG. 8B, since NLF remains substan- 
tially equal to zero from data points zero through almost 



300, the signal processor remains in the tracking mode. 
Further, since at no point prior to approximately data 
point 700 does NLF reach a value of 20, the signal proc- 
essor remains in the tracking mode until approximately 

s data point 700. 

During this long time period in which the methods 
and apparatus of the present invention maintain the sig- 
nal processor in the tracking mode, a variety of interfer- 
ing events are cause to occur in order to demonstrate 

10 the adaptability of the present invention. At time periods 
corresponding to data points surrounding data point 
100, liquid in tank 10 is agitated, as can be seen by the 
reduction in SNR at this point in FIG. 8A. However, uti- 
lizing the methods of the present invention of controlling 

is the masking thresholds and the window size, the target 
is never lost for more than a couple of spectral analysis 
repetitions, as can be seen in FIG. 8B. Throughout this 
agitation, the signal processor provides a substantially 
constant range indication of approximately 2.5 meters 

20 to the surface of the liquid in the tank. 

Beginning at the time period indicated by approxi- 
mately data point 200, the level of the liquid in the tank 
is increased by about 1 meter such that the range 
decreases from 2.5 meters to about 1 .5 meters. The 

25 change in the liquid level in the tank causes the SNR to 
decrease during this time period. With the increased liq- 
uid level causing a decreased SNR, the liquid in the tank 
was once again agitated at time period near data points 
300. As can be seen in FIG. 8B ( the combination of agi- 

30 tation with a reduced SNR caused by high liquid levels 
results in NLF taking on values greater than zero, thus 
representing a temporary loss of the target during one 
or more repetitions of the spectral analysis. However, 
since the methods of the present invention allowed the 

35 target to be quickly found, signal processor 110 does 
not leave the tracking mode as can be seen in FIG. 8D. 
Further, during this entire time period, the level of the 
liquid in tank 10 is accurately measured as is indicated 
in FIG. 8C. 

40 Next, beginning at a time period corresponding 
approximately to data point 350. the level of the liquid in 
tank 10 is once again decreased by about one meter. 
During the time period corresponding to between about 
data point 350 and about data point 600, water was 

45 sprayed into antenna 26 to provide interference. At the 
same time, during the time period between about data 
point 500 and about data point 600, the liquid in the tank 
was also agitated. Both of these interfering events are 
illustrated in FIG. 8A. While spraying water into the 

so antenna alone did not result in even a temporary loss of 
the target, the combination of spraying water into the 
antenna while agitating the liquid caused the target to 
be lost during tracking mode for a few repetitions, as 
illustrated in the NLF plot of FIG. 8B. However, during 

55 this entire time period, signal processor 110 remained 
in tracking mode as illustrated in FIG. 8D and accurately 
determined the level of the liquid in the tank as illus- 
trated in FIG. 8C. At a time corresponding to about data 
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point 600, the level of the liquid in tank 10 was once 
again increased. 

Beginning at the time period corresponding to 
about data point 600, the level of the liquid in the tank 
was once again increased by about 1 meter, thereby 
again dropping the SNR due to the corresponding 
increase in near field noise, as can be seen in FIG. 8A. 
Between the time periods corresponding to about data 
point 700 and about data point 800, the peanut oil in the 
tank was agitated while water was simultaneously 
sprayed into antenna 26. As can be seen in FIGS. 8B 
AND 8D, the combination of a very high liquid level, a 
low dielectric material (peanut oil), agitation and spray- 
ing water into the antenna caused the target to be lost 
for more than 20 spectrum analysis cycles. Thus, signal 
processor 110 entered the search mode of operation 
and remained there until a target was found for eight 
consecutive spectrum analysis cycles. 

Claims 
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20 



1 . A method of operating a FMCW radar tank level 
gauge that measures a level in a tank by obtaining 
a set of data points of mixed transmitted waves and 
received waves, the method comprising: 2 s 



generating and storing a set of spectral data 
values including a received target marker indi- 
cating the level; 

calculating an adaptive ^et of masking thresh- 
old values corresponding to at least a portion of 
the spectral data values; 
comparing the adaptive set of masking thresh- 
old values with the corresponding spectral data 
values to identify at least one spectral data 
value associated with the level; 
calculating the level of the tank from said at 
least one spectral data value; and 
reporting the level of the tank. 



30 



35 



2. 



3. 



The method of claim 1 wherein the step of calculat- 
ing the adaptive set of masking threshold values 
includes adapting individual masking threshold val- 
ues as a function of the associated spectral data 
value and at least one adjacent spectral data value. 

The method of claim 1 or 2 wherein the step of cal- 
culating the adaptive set of masking threshold val- 
ues includes adapting individual masking threshold 
values as a function of the associated spectral data 
value and past masking threshold values. 



40 



45 



5. 



searching for an absence of a later target 
marker in the window; and 

increasing the window to a second larger band- 
width if the target marker is absent from the 
window. 

A FMCW radar tank level gauge that measures a 
level in a tank by obtaining a set of data points of 
mixed transmitted waves and received waves, par- 
ticularly for carrying out the method according to 
any one of the preceding claims, the gauge com- 
prising: 

means for generating and storing a set of spec- 
tral data values including a received target 
marker indicating the level; 
means for calculating an adaptive set of mask- 
ing threshold values corresponding to at least a 
portion of the spectral data values; 
means for comparing the adaptive set of mask- 
ing threshold values with the corresponding 
spectral data values to identify at least one 
spectral data value associated with the level; 
means for calculating the level of the tank from 
said at least one spectral data value; and 
means for reporting the level of the tank. 



50 



4. The method of any one of the preceding claims and 
further comprising: 

defining a window of a first bandwidth around a 
received target marker as a function of the 
comparing; and 
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